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Abstract 

We study the asymptotic behavior of nonnegative solutions of the semilinear parabolic 
problem 

ut = Au + uP, X G M^, t>0 
u{0) = no, X G M^, t = 0. 
It is known that the nonnegative solution u{t) of this problem blows up in finite time 
for 1 < p < 1 + 2/N. Moreover, if p > 1 + 2/N and the norm of uq is small enough, 
the problem admits global solution. In this work, we use the entropy method to 
obtain the decay rate of the global solution u{t). 
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1 Introduction 



We considerer the semilinear parabolic problem 

ut = Au + uP, xeR^, t > 

(1) 

m(0) = Mo, xe M^, t = 0, 
where p > 1 and uq is nonnegative and nontrivial. The interest in this article is to study 
the asymptotic behavior of global in time solutions in order to obtain their decay rate. 
Although this matter has been treated for other authors, we do it using the entropy 
method. This method has been successfully applied, for instead, by J. A. Carrillo and G. 
Toscani [Jj for the asymptotic behavior of global solutions of some Fokker-Planck type 
equations. 

In [3], Fujita considered the evolution problem ([1]) and proved the existence of a critical 
exponent p* = 1 + which is called the Fujita's exponent. This exponent satisfies 

• for p > p*, if the norm of Uq is small enough, there exists a classical global in time 
positive solution of ([1]) which decays to zero when t ^ oo, in the other case the 
solution blows up in finite time; 

• for 1 < p < p* and any choice of the initial condition uq, every positive solution of 
([1]) blows up in finite time. 

In the case p = p* Hayakawa [5] proved that every positive solution of ([1]) has the same 
behavior as 1 < p < p*. This fact was proved by Kobayashi, Sirao and Tanaka [8] too. 

The aim of this work is to obtain the decay rate of the global solution u{s) of the 
problem ([T]) when t ^ oo. More specifically, we deduce that under certain hypotheses, 
the L^-distance of u{s) decays with the rate 

N 

\Ht)\\L2 < c {t + 1)'- , t>h, 

for certain time ti > 0, as well as the L^-distance of u{s) 

II / X II 1 2(N 1 \ 

The organization of this work is the following. The section 2 contains the notation 
and some previous results which are necessary for the next sections. In section 3 a 
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short summary about the entropy method steps is given. In the section 4 we obtain the 
exponential decay for the entropy production and for the entropy. The last section is 
devoted to deduce the asymptotic behavior of the mentioned global solution. 



2 Preliminaries 

The idea is to transform the equation in problem ([1]) in order to obtain significant in- 
formation on the asymptotic behavior of the global solutions. The change of variables 
utilized for Kavian [6] and Kawanago [7] is here employed. That is, 

x = {t+ ly/^y and t = e' - 1. (2) 

Then, v{y, s) results to be a solution of the problem 

Vs = Av + ^.Vv + ^ + vP, yeR^, s>0 

* 2 p—1 ' T ' ^g-j 

Let us observe that problem ([3]) has the same initial condition as problem ([1]). Moreover, 
we work with the following spaces 



Hi = {fe Lll V/ G L]] , Hi = {fe Hi/ V/ G Hi) , 
where p{y) = exp(|yp/4) and r > 1 is a constant. Related to these spaces are 

(/,^?)l?= / fgpdy, ||/|U. = (/,/)^ 

{f,g)Hi = {f.g)Li + {Vf,Vg)Li, 
ll/llL + l|V/|li.; 



HI 



Observe that the equation in ([3]) has to be written in the shape 

Vs = -Lv H + f ^ , 

P — 1 
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where L is the self-adjoint operator given by 



Lv = —Ay — 2'^^' defined over D{L) := Hp. 



We know that this operator satisfies 



Ai = Y is the least eigenvalue of L. Then, the following inequality holds 



N „ „ 
2 " " " - 



|Vf l^p dy 



(4) 



the operator L has compact inverse. 



As well, it is known that u^o = C{N,p) \x\ p-i with C{N,p) 
singular equilibrium of ([1]), that is, Uoo is a solution of the Lane-Emden equation 



p-i \ p-i 



IS a 



Au + uP = 0, xeR^ u>0, N>3 



which arises in astrophysics and Riemannian geometry. It is well-known that this fact is 
possible only for those values of p that verify p > j^z^, since Gidas and Spruck [1] proved 
that there are not stationary solutions in any other cases. In 1993, Wang [9j proved that 



if > 3, p > 



N 
N-2 



and 



< Uo{x) < A Moo, 

where < A < 1, then ([T]) has a unique global classical solution u with < n < Xuoo- It 
also satisfies that 

u{x,t) < [{X^-P -l){p-l)t]'^K 
This inequality, in terms of the problem ([3]), can be expressed as 



1 



e'' - 1 



(5) 



[(Ai-P-i)(p_i)]^ 

/^From now on, we assume that the Wang's theorem hypotheses are satisfied. Moreover 
the framework will be the set of global solutions f of ([3]) such that v & X, with 

X=(/Gi/;nL-/ />Oand lim/ \V f {s)\^ p dy = o] . 
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3 Entropy method 

We study the asymptotic behavior of the global in time solutions of the problem ([3]). 
For it, we use the already mentioned entropy method. The essential application of this 
method will consist in the following steps. 

• Define a suitable entropy functional E{v{s)) for the equation ([3]) and study its 
properties. 

• Compute the entropy production 

I{v{s)) = ^^E{v{s)) . 

• Compute the derivative of entropy production and obtain a differential equation of 
type 

^^I{v{s)) = -CI{v{s))~R{s), 
for certain constant C > and certain function R[s). 

• Prove the properties of R{s) that permit to obtain an exponential decay of I{v{s)), 

I{v{s)) < A e-^' . 

• Obtain the same decay rate for E{v{s)) from the previous items, more specifically 

E{v{s)) < B e~^'^ , for s > Si and certain si > 0. 

• Give a bound of ||f ||l2 in terms of the entropy and entropy production which permits 
to get conclusions on the decay of the mentioned norm. 

The same entropy functional introduced by Kavian and Kawanago (|6] and [7] respec- 
tively) will be used in the present article. 

Definition 3.1. For every v E fl L^"*"^ the entropy functional is defined by 

E(v) = I I \Vv\' - ^. ^ ^. - v^^'] p dy . 

2{p-l) p+l J 
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In order to obtain the decays announced above, some properties of this functional are 
needed. These properties are summarized in the next proposition. The first of them was 
proved in [B]. 

Proposition 3.2. Let uq G Hj^ fl L°°,Uo > 0, E{uq) < oo and v = v{y,s) the global 
solution of f G X. Then 

1. if there exists sq > such that E{v{so)) < and v{so) 0, v blows up in finite 
time; 

2- iE{v{s)) = -I{v{s)) where I{v{s)) = j^^v^p dy; 

3. there exists M := hm E{v{s)) and, moreover, M = 0. 

Definition 3.3. The functional I {v{s)) of the proposition \3.2\ is called entropy produc- 
tion. 

Proof. For the first property, the reader is refered to [6]. The second one is deduced by 
derivating E{v{s)) with respect to s, integrating by parts and keeping in mind that v 
is solution of ([3]). In order to prove the third property, first observe that owing to the 
second one ^we have that E{v{s)) is non increasing. As < E{v{s)) < E{v{0)) < oo, the 
existence of the limit is warranted. Therefore, we must only see that the limit is equal to 
zero. For this issue, we observe that the first term in the expression of E{v{s)) tends to 
zero when s — > cxd since v E X. The second term of that expression goes to too when 
s oo due to the inequality (jlj). To see the behavior of the third term of the entropy 
functional we use the inequality ([5]) and obtain 

1 



< / ^—v'P+^pdy 

Irn p+l 



< 



-1 



1 



(Ai-p-l)(p2-i) 



v'^p dy. 



The last expression can be bounded for large values of s as following 

^ -v^+'p dy < C-—-4rrT-^ [ ^'P dy, 



p+1 ^ ^ - (Ai-P-l)(p2-l^ 
where C is a positive constant. Taking into account the last inequality, we get that the 
third term of E{v{s)) tends to zero when s — > oo. □ 



Now, we want to prove the decay of the entropy production and, as a result, the decay 
of the entropy functional. It will be made in the next section. 



4 Decay of the entropy functional 

The computation of '^^^'^'f^^ is needed to obtain the decay of the entropy production. We 
write this derivative in a convenient way using that v is the solution of the equation ([3]). 
That is, 



—I{v{s)) = / 2vsVsspdy 
as J^N 



^ livis)) - 2 (Lv,, V,) + 2p I v'-\Ip dy 



p — 1 

= -2jI{v{s))-2R{s), (6) 

where 7 is a positive constant and R{s) is an appropriate function. Both 7 and R{s) are 
defined by 

N 1 

7 = V 7>0, 

2 p — 1 

N f 

R{s) = {Lvs,Vs)lI \\vs\\l2 -p v^-^vlpdy. (7) 

We need to have more information about R{s) that allows to obtain conclusions about 
the decay of the entropy production. This is the aim of the next three lemmas. 

Lemma 4.1. Let v G and Q C an open set. Then 

f N 

/ {vLv — — t>^)p dy > . 
Jn 2 

Proof. We observe that, if Xq is the characteristic function of the set Q, then 

/ vLv p dy = / (vLv) Xq p dy 

V Xn{Lv) Xq p dy . 
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Owing to (Lv) Xq = L{v Xq) almost everywhere, we have that 

/ vLv p dy = / V XqL{v Xq) p dy 



N 

> - I {vXnfpdy 



v"^ Xq p dy 



2 

N 

y 

= y v pdy, 

where the inequahty in ([8]) is a direct consequence of (jlj). From the last computation the 
lemma's statement is established. □ 

Lemma 4.2. Let f E L^(R^) such that J^f dy > for every open set C M^. Then, 
/ > almost everywhere in M^. 

Proof. We suppose that / isn't nonnegative a.e., then there exists a measurable set Q 
such that m{Q) > and / < in fi. Here, m{Q) denotes the Lebesgue measure in of 
the set Q. Then, 

a = / f dy < 0. 
Jq 

For each n eN, there exists an open set Gn C such that Q C Gn and m{Gn — Q) < K 
We can choose those open sets (?„ in such a way that the sequence {Gn)n is increasing 
with the inclusion. Then, we observe that 

f dy = a+ f dy 

Gn ^ Gn — f2 

= a+ ipndy (9) 
Jrn 

where ipn = X(^G„-n)f- It is quite clear that the functions ipn are integrable functions, 
they satisfy \(pn\ < |/| and, moreover, they verify that when n oo a.e. Then, 

owing to dominated convergence theorem, it results that 

lim / ifindy = 0. 

Thanks to the last equality we deduce that there exists a natural Nq, such that for every 
n>No 
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< 



\a\ 



(10) 



Using this inequality in ([9]), we see that it verifies 



f dy < a + 



Gn 



I 



dy 



\a\ a 



for every n > Nq. That is, the integral of the function / on the open sets Gn, for n > No, 
is negative. This fact contradicts the lemma's hypothesis. □ 

The next lemma provides a bound of R{s) in terms of a new nonnegative function 
with exponential decay. 

Lemma 4.3. Under the same hypotheses in the Wang's theorem and if, moreover, A < 



3p-l 



or 



^(p_T)2 + 1 J then the function R{s) defined by ([^ satisfies R{s) > —\ K{s) f 
certain function K{s) that verifies 

1. K{s) > 0, 

2. there exists a constant a > and a time Si > 0, which depend on p, N and X, such 
that 

K{s) < K{si) e-(27+«)^^ for s > si, and 

3. Jq°° e'^'^^K{s) ds < C , for a suitable positive constant C . 
Proof. We define the function K{s) as follows 

K{s) = 2p I v^-^vIp dy. 

It is clear that K{s) is nonnegative. To see 2, we compute the derivative of the function 
K{s) and write it conveniently. 

dK{s) 



ds 



2p I [{p - 1) yP-^l + 2 vP-'vsVss]p dy 

^ K{s) + 2p {3p - 1) / v^'^'\'i p dy 



p — 1 



—2p (p — 1) / '^VgLv p dy — Ap I ^VgLvs p dy. 



(11) 



To get a bound of the second of the four terms of (fTTI) . we use the inequahty ([5]) in order 
to obtain that 

2p (3p - 1) I v^'-V. p ,y < J! ^ K(s). (12) 

For the bound of the third term, we apply first the lemmas [4 . 1 1 and in order to deduce 
that vLv — Y f ^ > a.e., recalling that v'''~'^Vg > we conclude 

{vLv - Y'i^^) v^'^v^ p dy >0 . 

This inequality quickly leads to a bound of the third term of ffTTj) which is a multiple of 
K{s). 

2p (p - 1) / yP-^v'^^Lv pdy > 2p (p - 1) ^ [ v^-^v^, p dy 

= {p-l)^K{s). (13) 

For the last term of ( ITTi) . we use first that VgLvg — y "wf > a.e. (it's owing to the lemmas 
Oandiil]) and that also v'P'^ > 0. We obtain 



{vsLvs-vD vP~^ pdy>0. 

This inequality permits to bound the last term of f|TT]) by a multiple of K{s) as follows 

Ap / v^'^VsLvs p dy > 2p N v^'^vl p dy 

= N K{s). (14) 



Then, taking into account f[T2l) . f[T3|l and flMj) . we have that 

K{s) 



ds \ p — 1 

3p — 1 e' 



where 



(Ai-P - l){p- 1) - 1 
-27 K{s)~p{s) K{s). 



p{s) = -l + {p-l)- ^ ^ ' 



2 {X^-P- l)(p- 1) -1 
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with /(s) = ^FTT- Now, we must prove that > for s > si, where si is a time which 
depends on p, N and A. It is clear that, owing to (fT5l) . it is equivalent to require that for 
/(s) < B, where B is the constant defined by 



B 



3p-l 



We can take Si = f ^(^-^) in order to get > for s > si, since -B > 1 due to the 
hypotheses of the lemma. Then, a bound of the derivative of K{s) is obtained from 

dK{s) 



ds 



<-{2-f + a)K{s), s>si. 



where a is a positive constant (we can take, for example, a = ^^^). Thus, from a certain 
si > 0, the function K{s) has exponential decay rate. That is, 

K{s) < K{si) e-(2^+'^)^ (16) 

This statement proves the second part of the lemma. The last part of the lemma can be 
deduced immediately from the previous one. □ 

Applying the former lemma to ([6]), we obtain the mentioned decay for the entropy 
production as it will be proved in the next theorem. 

Theorem 4.4. Under the hypotheses of the lemma \4^ and if, moreover, I{uq) < oo, then 
I{v{s)) has an exponential decay rate. More precisely, 

Hvis)) < [Iiuo) + C] e-2T^ 

Proof. From the expression ([6]) and the lemma 14. 3[ we have 

{e''^' I{v{s))) = -2 e^^' R{s) < e^^' K{s). 

Integrating between and s we obtain 

e^^' I{v{s))-I{v{0)) < [ e^^^K{a)da 

Jo 

POO 

< / e^^" K{a) da<C 
Jo 

.: I{v{s)) < [I{uo)+C] e-2^^ 



□ 
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Now, we are in a position to prove that the entropy functional decays exponentially. 
It will be proved in the next theorem. 



Theorem 4.5. Under the hypotheses of the theorem 4-4- the entropy functional E{v{s)) 
has exponential decay, that is, 

E(v{s)) < C e~'^^^, for every s > si, 

for certain Si > which depends on p, N and X, and for certain positive constant C that 
depends on I{uq), p, N and A. 

Proof. From part 2 of proposition 13.21 and expression IQ, we have that 

dEivjs)) ^ 1 dlivjs)) ^ 1 ^ 
ds 27 ds 7 

> ±^M^-±K(s), (17) 

- 27 rfs 27 ^ ^' ^ ^ 

where the inequality in (fT71) is owing to the bound of R{s) from K{s). Now, we integrate 
between s and b in f|T7j) and we use part 2 of lemma 14.31 to obtain 



E{v{b)) - E{v{s)) > 



1 r^/_v,_^^ 1 f K{si) _(2^+„y 



> TT- [livib]) - I{v{s))] - - -^e 
Z'-f Z'-f \27 + a 



b 



for every s > Si. Taking the limit for 6 — > oo and using part 3 of proposition 13.21 we get 

E(v(s)) < ^ + — P^i^ e-(27+-)« 
^ ^ - 27 ^ ^ 27 (27 + a) 

for every s > Si. Therefore, owing to the theorem fl4.4l) . the announced decay for the 

entropy functional takes place 

E{v{s)) < C e"^''"', for every s > Si. 

□ 



5 Asymptotic behavior of the solution 

To finish with the application of this method, we must obtain a bound of the norm of the 
solution v{s) in terms of the entropy E{v{s)) and the entropy production I{v{s)). For 
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this purpose we define the following function g{s) 

9{s) = \ I v'^{s) p dy . 



2 

The next lemma provides a bound of the function g{s) in terms of E{v{s)) and I{v{s)), 
that is, a bound of the norm of v{s) in the space L^. 

Theorem 5.1. Under the hypotheses of theorem o,''T'd if, moreover, p > p , where p 
is defined by 



N 



^ = 3, 



1 + 1, N>A 



p = 
then 

1, [S^N- 2) g{s) < \ I{v{s)) + (p + 1) E{v{s)) and 

2- g{s) < C e~'^"'^ for all s > Si, where Si is a positive number which depends on p, 
N and \, and C is certain positive constant which depends on I{uo), p, N and A. 

Proof. We observe that the derivative of the function g{s) satisfies 

g\s) = / vvs p dy (19) 



< 



1 2^1 2 

- V + - v: 



pdy = g{s) + ^I{v{s)). (20) 



2 2 

On the other hand, we can obtain another expression for g'{s) replacing Vg in (fT9|) accord- 
ing to the problem ([3]) and using the definitions of E{v{s)), L and g{s), that is, 

g'is) = -ip+l) E{vis)) + ^ {Lv, v) - g{s). (21) 

Owing to dH), (I20D and ([21]), we have that 

N-2] g{s) < I I{v{s)) + {p+l) Eiv{s)) . 



2 y ' - 2 

As p > + 1 and using the bounds for the entropy and entropy production in theorems 
14.41 and 1^31 respectively, we get 

g{s) < C e"^'*'* para s > si , 

where C is a positive constant which depends on A^, p, A and /(f (0)). □ 
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Thanks to the definition (fTSi) . we can deduce the decay of the norm of v{s) in the 
space which is the aim in the following theorem. 

Theorem 5.2. Under the hypotheses of the theorem \5.1[ then 

\Hs)\\Li<Ce-^\ s>s,, 

for si > which depends on p, N and X, and for certain positive constant C which 
depends on I{uq), p, N and A. 

Now, we can use again the change of variables ([2]) in order to obtain the decay of the 
solution u{x,t) of ([I]), getting the next result. 

Theorem 5.3. Under the hypotheses of the theorem \5.1\ then 

N 

Mt)\\L2 <c{t + i)-- , t>ti, 

forti > which depends onp, N andX, and for certain positive constant C which depends 
on I{uq), p, N and A. 

Remark 5.4. Let us observe that, because of the Wang's theorem and the theorem \5. 31 
we obtain the decay of the norm \\u{t)\\Lq for q > 2. This is, 

1 2 { N 1 \ 

\Ht)\\Ll r^t'~'~^^~'~K (22) 

for t>ti. 

Remark 5.5. Notice that, although the decay of the norm \\u\\li ~ t ^ I q) obtained by 
Kawanago in /?/ is better than [2^) in the case q > 2 and the same for q = 2, the result 
obtained in ^ is true for < p < and the result obtained in this work corresponds 
to the range p > p . 
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